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Abstract 



We present an overview of the phenomenology of prompt neutrinos produced in the 
atmosphere by the decay of charmed particles. In face of the great activity involving 
operational and proposed large scale neutrino telescopes, we discuss the implications 
of different scenarios in calculating the prompt flux, upon the detectable signals in 
these experiments. 
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1 Introduction 



One could imagine to build a special kind of experimental facility to study 
high-energy and very-high-energy neutrinos. While on the wishful thinking, 
one could gather all the best conditions. Desirable characteristics include: 

• to provide acceleration of parent particles up to the highest energies; 

• to be coupled to the largest particle beam dump; 

• to be shielded by a set of the largest and most dense absorbers available; 

• to be monitored by an ensemble of several detectors, each operated with a 
different detection technique, and complementary acceptances; and 

• to produce and detect neutrinos covering a broad range of energies, perhaps 
from lO^eV to lO^^eV. 

Does it sounds like a too difficult enterprise? Well, we get very close to it, 
if we regard all cosmic-ray phenomena and corresponding experiments in a 
complementary way. 

Think of several cosmic particle acceleration mechanisms, provided by sources 
such as supernova explosions, neutron stars in binary systems, active galactic 
nuclei, gamma-ray bursts, pulsars, rapidly spinning magnetars and yet some 
exotic objects (topological defects, decaying massive relics). Choose as target 
our planet, where gravitationally trapped gaseous elements serve as a beam 
dump, and the solid core material, from one surface to the other, provide at 
the same time shielding and support for different purpose detectors, scattered 
around at different heights and depths. A block diagram of this "ultimate" 
high-energy facility looks like: 
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Fig. 1. Sketch of the "ultimate" experimental neutrino facility (the Cosmos in- 
cluded!). 

A more realistic sketch for this omnipresent experimental neutrino facility 
is presented in Figure 1. One should consider the whole cosmos as part of 
the figure. One may also include several detectors scattered around, with 
different and complementary purposes, for example: optical and gamma-ray 
space telescopes; primary cosmic-ray satellite and balloon borne calorimeters; 
mountain-altitude cosmic-ray family emulsion chamber detectors; cherenkov, 
fluorescence and particle extensive air shower detectors; muon magnetic spec- 
trometers; and yet large underground, underwater and under-ice observatories. 

High-energy extraterrestrial neutrinos (above tens of TeV) may be directly 

produced at the same (galactic or extra-galactic) sources of high-energy cosmic- 
ray particles and gamma-rays, including those of the highest energies. Low en- 
ergy neutrinos (below tens of MeV) are produced directly in nuclear reactions 
at the core of the Sun. 



3 



Intermediate energy neutrinos are produced locally, at the planetary scale. A 
cosmic accelerated primary particle, a miclcon for example, impinging upon 
nuclei of air in the atmosphere, will generate a shower of secondary particles, 
which may propagate to the ground or even below. The main contributions 
to air showers come from the prodution of a) neutral pions, giving rise to 
electromagnetic cascades (a sequence of photons and electron-positron pairs); 
and of b) charged mesons (mainly pions and kaons), giving rise to hadronic 
cascades. Charged-meson decays will give rise to neutrinos and muons; these 
muons on their turn may eventually decay into neutrinos, all contributing to 
what is called the conventional atmospheric neutrino fluoi{l]. 

The semileptonic decay of heavy quark particles, produced in the hadronic 
cascades initiated by cosmic-ray interactions, also gives rise to atmospheric 
neutrinos, a particular set known as the prompt atmospheric neutrino flux. 

These prompt neutrinos are the subject of the present talk. Among the moti- 
vations for their study, we highlight the following: 

• High Energy Neutrino Telescopes - 

Though neutrinos are the most abundant cosmic rays at sea level, their 
measurement requires large underground detectors, in order to compensate 
for the small neutrino-nucleon cross section. Nowadays operational neu- 
trino telescopes, like AMANDA [2] and BAIKAL [3], arc taking and analizing 
data. They are being optimized to increase signal, data statistics and en- 
ergy coverage, and to reduce the background, mainly caused by muons and 
by electronics limitations. Other similar and much larger instruments are 
being contemplated (Antares, Nestor, Nemo, lceCube)[4]. As will be shown 
later, there is the possibility that prompt neutrinos are observable in such 
neutrino detectors [5]. 

• Background to Cosmic Neutrinos - 

The measurement of galactic and extra-galactic neutrinos, which is the main 
purpose of neutrino telescopes, is of course limited by discrimination against 
atmospheric neutrinos, the prompt component included. 

• Bounds to Charm Production Cross Section - 

The actual measurement of prompt neutrinos, or the failure to do so, will 
impose constrains to the charm production cross section at high energy [6], 
in a similar way as it has been done through the investigation of muon poor 
horizontal air showers[7]. 

• Proton Small-x Gluon PDF - 

It has been shown that the spectral index of the prompt atmospheric neu- 
trino flux may depend linearly on the slope of the gluon distribution function 
at very small x, not reachable at colliders [8]. 

• Background to High Energy Neutrino Oscillations - 

Prompt decay of the charm-strange meson Dg provide the unique direct 
source of high-energy tau-neutrinos in the atmosphere, which may compete 
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with those coming from flavor oscillation of atmospheric or extraterrestrial 
muon- neutrinos [5] . 

Why charm turns out to be of such particular interest to the study of high 
energy neutrinos? 

The answer comes easily from the analysis of the critical energy of different 
neutrino-parent particles produced in cosmic-ray interactions in the atmo- 
sphere. The critical energy delimits the competition between decay and inter- 
action lengths, that is to say, above this energy the parent particle is likely 
to interact or be slowed down, rather than to decay into a neutrino. Table 1 
compares the critical energies calculated for muons, pions and kaons (parents 
of the conventional atmospheric neutrino flux) and for relevant charm baryons 
and mesons (parents of the prompt flux) . While the contributions of the con- 
ventional component to the neutrino flux are constrained below 1 TeV, the 
charm prompt decays - although less copious - represent the only atmospheric 
contribution up to around 100,000 TeV. 

In what follows we present an outline of the calculation of the prompt neu- 
trino component in the atmosphere, comparing different phenomenological 
approaches found in the literature. Next wc discuss the implications of the 
different prompt flux scenarios upon the detectable signals in large scale neu- 
trino telescopes. 



Table 1 

Critical energies ScriUc of selected particles. 



Particle 


Elementary 




^critic 




contents 


(MeV) 


(GeV) 


^+,^- 


lepton 


106 


1.0 


7r+, 7r~ 


ud, ud 


140 


115 


K+,K- 


us, us 


494 


855 


A" 


uds 


1116 


9.0 X 10^ 


D+,D- 


cd, cd 


1870 


3.8 X lO'^ 




cu, cu 


1865 


9.6 X 10'^ 




OS, cs 


1969 


8.5 X 10^ 


A+ 


udc 


2285 


2.4 X 10^ 
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2 Calculation 



The particle and energy flow along the development of an air shower, at differ- 
ent atmospheric depths x (given in units g/cm^), relevant for the production 
of prompt leptons, can be summarized as: 



Primary 
Cosmic Ray 



-(7+1) 
N 



Nucleonic 
Flux 



Charm Production 
and Decay 



d^i/dx {Ei,x) 



i 



Prompt Lepton 
Flux 



HEi,xi) = Jdx JdE, ^ B,j:<^>,{E,, 



X) 



The index stands for the nucleons (protons, neutrons) arriving at the top 
of the atmosphere {x = 0), with differential energy spectrum $iv described 
by a power law. They cascade down, subjected to an attenuation length A^, 
to a generic level x where they may produce a charmed particle of type-i 
(i = D^, D°, D°, Df,A^). The term Ki is related to the charm production 
spectrum- weighted moment Zj\fi, while Aj and di are respectively the charm 
interaction and decay lengths. Finally, upon decay of the charmed particle to 
produce the lepton I, with branching ratio and decay spectrum df/dEi, we 
have calculated the differential flux of prompt leptons ^i{Ei,xi), with energy 
El at atmospheric depth xi. 

Detailed calculation of the prompt lepton component is presented in Ref. [9]. 
In order to have a taste of what this result may look like, and what it signifles, 
we write down an approximate solution, valid for energies below the charm 
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critical energy and for depths of the the order of the whole atmosphere: 

Prompt Lepton <— Decay <— Produce <— Attenuate <— Primary Flux 

Reading from right to left, this expression tells us what we have seen in the 
previous diagram: the primary flux is attenuated, then produces charm parti- 
cles which decay into prompt leptons (the decay spectrum-weighted moment 
Zii is related to the decay spectrum df /dEi). 

Does it seems that simple? Well, we must first take a look at the many in- 
gredients used in the calculation, considering that there are many possible 
phenomenological assumptions to be adopted for each. 

3 Phenomenology 

A preliminary list of parameters involved in the cascade routines leading to 
the calculation of the prompt flux comprises: 

• Primary Spectrum Amplitude No and Slope 7 

• Nucleonic Inelastic Cross Section of^"''^^(E) 

• Nucleonic Interaction Length Xn{E) 

• Nucleonic Spectrum- weighted Moment ZNN{'y,E) 

• Nucleonic Attenuation Length An{E) 

• Charm Inelastic Cross Section (/^^^"'""(E^) 

• Charm Interaction Length \i{E) 

• Charm Production Spectrum- weighted Moment ZNi{'y,E) 

• Three-body Semileptonic Decay Spectrum ^ 

• Charm Decay Spectrum- weighted Moment ^^(t) 

• Charm Decay Length di 

It is up to the user to choose a desired brand for each ingredient. As an 
example, consider some charm production models found in the literature: 

Quark Gluon String Model (QGSM) - A semi-empirical model of charm 
production based on the non-perturbative QCD calculation by Kaidalov and 
Piskunova[10], normalized to accelerator data, and applied to the prompt 
muon calculation by Volkova et a/. [11]. 

Recombination Queirk Parton Model (RQPM) - A phenomenological 
non-perturbative approach, taking into account the contribution of the in- 
trinsic charm to the production process, in which a cc pair is coupled to 
more than one constituent of the projectile hadron, as described by Bugaev 
et a/. [12]. 
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Perturbative QCD (pQCD) - There are several calculations based on per- 
turbative QCD, for example Thunman et a/. [13] evaluate explicitly the 
charm production up to leading order (LO) in the coupling constant, us- 
ing the Monte Carlo program PYTHIA, and include the next-to-leading 
order (NLO) distribution effects as an overall factor. More recently Gelmini 
et a/. [14] updated the calculation to include the full contribution of NLO 
predictions to the Icpton fluxes. Pasquali et a/. [15] offer an alternative im- 
plementation of the perturbative QCD approach. 

Limiting Conditions - E. Zas et a/. [16] calculate extreme cases of charm 
production, at both low and high production rate limits. Higher rates are 
achieved assuming a charm production cross section which is 10% of the 
total inelastic cross section (called Model- A), extrapolated with a log^(s) 
energy dependence. Lower rates are obtained using pQCD at NLO, with 
structure functions given by Kwiecinski-Martin-Roberts-Stirling and adopt- 
ing relatively hard parton distribution functions (called Model-E). 

Prompt muon fluxes calculated assuming these charm production models are 
compared in Figure 2. (The fluxes are multiplied by in order to flatten the 
curves.) Two calculations of the conventional atmospheric muon flux are also 
included, and we note that the cross-over of the conventional and prompt com- 
ponents may occur at any energy between 2 x 10^ GeV and 2 x 10^ GeV. The 
prompt flux intensity, itself, spans up to four orders of magnitude! However, 
we call attention to the fact that the comparison is somewhat unfair, since 
not only different charm production models are being compared: calculations 
of different authors involve different assumptions for the other cascading in- 
gredients, so that we are seeing a combination of effects. 

Ref . [9] discusses the influence of each separate ingredient upon the final prompt 
fiux, adopting diverse parametrizations available in the literature. For exam- 
ple, the solely effect of the primary cosmic-ray spectrum is displayed in Fig- 
ure 3. Curves in a) show the spectrum at the top of the atmosphere, according 
to 5 models. Curves in b) are the corresponding prompt muon-neutrino ver- 
tical fiux at sea level. Same comparison can be made, fixing all ingredients 
but the charm production spectrum- weighted moment. The resulting muon- 
neutrino fiuxes for 5 different parametrizations of charm model are shown in 
Figure 4. 

We can think now of mixing different parametrizations of all ingredients listed 
previously. It is possible to choose extreme combinations, that is to say, those 
leading to lower and higher prompt lepton fluxes, for a given charm model. 
This procedure can alter the flux which has been reported in the literature: 
Figure 5 shows how the calculation of Volkova et al. can be shifted down and 
how the calculation of Thunman et al. can be shifted up, by adopting different 
ingredients, other then the charm model. With the same procedure. Figure 6 
shows how the calculation of Bugaev et al. wiggles! 
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Fig. 2. Comparison of calculated differential vertical atmospheric muon fluxes at 
sea level, as reported by several authors (see text). 
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Fig. 3. a) Comparison of primary cosmic-ray energy spectra, as given by different 
parametrizations ; b)comparison of corresponding prompt neutrino fluxes, assuming 
all the ingredients fixed but the primary spectrum. 



Prompt Neutrino Flux with Charm Model: 

QGSM 

RQPM (+ Feynman Scaling) 

RQPM (+ ScaUng Violation) 

pQCD (no knee) 




lO' lO' lO' lO' lO' lO' lO' 
Neutrino Energy (GeV) 



Fig. 4. Comparison of prompt neutrino fluxes for different charm Z-Moment models, 
assuming all the other ingredients fixed. 
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Fig. 5. Prompt neutrino flux calculated with QGSM (upper band) and with pQCD 
(lower band). The bands represent the range between extreme ingredient combina- 
tions, for each charm model. Prompt flux from Volkova et al. (dotted line) and from 
Thunman et al. (dot-dashed line) also shown for illustration. 




Fig. 6. Prompt neutrino flux calculated with RQPM. The band represents the 
range between extreme ingredient combinations. Prompt flux from Bugaev et al. 
(dot-dashed line) also shown for illustration. 
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4 Implications 



We have seen how the prompt flux calculation is subjected to large uncer- 
tainties, resulting from the imprecise knowledge of both atmospheric particle 
showering parameters and the modeling of the production of charm. We ob- 
serve that the calculated neutrino fluxes may be shifted by up to one order of 
magnitude for a given charm model. When comparison is made among alter- 
native charm production models, the spread reaches two orders of magnitude. 
We will therefore define an allowed range between the maximum (MAX) and 
minimum (MIN) prompt neutrino fiuxes. The MAX combination of ingredients 
adopts the QGSM below the knee and the RQPM above, in addition to other 
extreme parameters. The MIN combination adopts the pQCD charm model, 
with ingredients that pulls the flux down (complete description of choices is 
listed in Ref.[9]). 

It is not difficult to extend the previous calculation to obtain the fiux of prompt 
tau-neutrinos, the main contribution coming from decay of D^-mesons into a 
rz/-pair, and subsequent in-chain tau decays producing more neutrinos[18]. 
The above mentioned conclusions on the uncertainties of the calculation still 
hold for the tau component. 

Applying the prescription of Gandhi et a/. [17], we can use the calculated fluxes 
(of prompt z/^ and z/^) to estimate the rates of neutrino-induced muons and 
taus, per year per effective detection area (in 27r sr), in operating and pro- 
posed neutrino telescopes. Table 2 summarizes the event rates for BAIKAL 
NT-200[3], which have an effective area of 2 x 10^ m^, for AMANDA-II[2], with 
effective area of 3x lO"' m^, and for two different thresholds in the proposed kni'^ 
experiment 1CECUBE[19]. The estimates for the conventional atmospheric 
neutrino component (CONV) are shown for comparison. Experimental sig- 
natures may come to the rescue on disentangling the two components. For 
example, the ~ 1,500 prompt z/^'s above 1 TeV in ICECUBE data may be ex- 
tracted using their flat zenith angular distribution (while the conventional flux 
has a strong angular dependence). Prompt z/^'s would produce characteristic 
showers in the detector, enabling discrimination of the ~ 200 (20) neutrino 
induced taus above 1 (10) TeV. 

The next question is how the conventional and prompt atmospheric fluxes 
compare to expected rates of extragalactic neutrinos[5]. Adopting for compar- 
ison the upper bound to the diffuse flux of high-energy neutrinos from cosmic 
accelerators, such as active galactic nuclei and gamma-ray bursts, we summa- 
rize the situation for muon-neutrinos in Figure 7. The bands are the allowed 
region for prompt neutrinos as in Figures 5 and 6. Also shown is the conven- 
tional atmospheric flux. Concerning extraterrestrial neutrinos, the two upper 
curves, labeled MPR (solid and dotted lines), represent the bounds imposed 
to neutrino fiuxes from sources that are optically thick or thin to the emis- 
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sion of neutrons[20], respectively. The two straight hnes labeled WB (dashed 
and long-dashed), represent bounds on sources which produce the highest en- 
ergy cosmic rays [21]. The higher flux allows for cosmological evolution of the 
sources. We observe, from Figure 7, that the prompt neutrino flux exceeds 
the conventional at some energy above 20 TeV. If the MPR bounds are to 
hold, the prompt component is merely a background to be discounted. How- 
ever, according to the WB-scenario, the prompt component may exceed the 
extragalactic diffuse flux up to 300 TeV, if the sources are to evolve, and up 
to 2 PeV if not. 

Allowing for flavor oscillation of atmospheric and extragalactic muon-neutrinos 
into tau-neutrinos. Figure 7 can be transformed into Figure 8, remarking that 
prompt Ut-'s are directly produced from Dg decays and are independent of os- 
cillation mechanisms. Once again we observe the possibility that there is an 
observational window opened for the prompt component, this time for tau- 
neutrinos, extending from 2 TeV up to 40 TeV (or up to 500 TeV, depending 
on the scenario). 

Table 2 

Upward-going niTion and tau event rates per year arising from charged current in- 
teractions of atmospheric neutrinos in ice or water, for different neutrino telescope's 
effective areas and thresholds. Prompt flux calculation based on maximum (MAX) 
and minimum (MIN) range limits described in the text, compared to conventional 
atmospheric flux (CONY). 







fi+ + n 




r+ 




Experiment 


CONV 


Prompt 


Prompt 


Prompt 


Prompt 






MAX 


MIN 


MAX 


MIN 


BAIKAL NT-200 


22 


3 











AMANDA-II 


330 


44 


2 


7 





ICECUBE (> 1 TeV) 


11000 


1470 


53 


216 


11 


ICECUBE (> 10 TeV) 


170 


157 


3 


22 


1 



5 Conclusion 

The calculation of the prompt lepton flux in the atmosphere is subjected to 
large uncertainties. Partially because of our imprecise knowledge of certain 
atmospheric particle showering parameters. Further uncertainty is associated 
with the extrapolation to high energy of a variety of models describing the 
accelerator data on charm production. Care must be taken when comparing 
results from different calculations, for the combination of ingredients may shift 
the fluxes up or down. 
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Fig. 7. Comparison of several contributions to the high-energy flux. The bands 
represent the allowed range between maximum and minimum atmospheric prompt 
neutrino fluxes for different charm production models. Thick dashed line is the con- 
ventional atmospheric flux. Solid, dotted, dashed and long-dashed lines correspond 
to upper bounds imposed to diffuse extragalactic neutrino flux by the observation 
of high-energy cosmic-ray and gamma-ray spectra, according to different scenarios 
explained in the text. Fluxes are multiplied by E'^. 



Once the uncertainties are quantified, in terms of an allowed band for the 
prompt atmospheric neutrino flux, we come to the result that prompt neutri- 
nos may be observable in a kilometer-scale neutrino telescope. Actually, it may 
be possible that a small component of prompt neutrinos are already present 
in the observed samples of the ongoing experiment AMANDA, given their 
current limit on the atmospheric neutrino flux[2]. Of course, the existence of 
a direct observational window for prompt neutrinos depends on the actual 
bounds to extragalactic neutrino fluxes and, in the case of tau-neutrinos, on 
the actual values of the flavor oscillation parameters, which are currently sub- 
jects of great interest and of intense research. 

The actual measurement of prompt neutrinos, and even the failure to do so, 
will certainly impose constrains to the charm production cross section, at 
energy domains not yet accessible to particle colliders [6]. 
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MPR-Thin 



Neutrino Energy (GeV) 



Fig. 8. Comparison of several possible contributions to the high-energy Vr flux. 
Again, the bands represent the allowed range between maximum and minimum at- 
mospheric prompt neutrino fluxes for different charm models. Thick dashed lines are 
the result of maximum mixing flavor oscillation of the conventional atmospheric 
flux, for Am^ around the Super-Kamiokande value 3.2 x 10~^ eV^ (thick dot-dashed 
line). Solid, dotted, dashed and long-dashed lines correspond to the upper bounds 
in Figure 7, subjected to vacuum flavor oscillation, averaged in transit to Earth. 
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